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Abstract 

We use covariant and first-order formalism techniques to study the properties of general 
relativistic cosmology in three dimensions. The covariant approach provides an irreducible 
decomposition of the relativistic equations, which allows for a mathematically compact and 
physically transparent description of the 3-dimensional spacetimes. Using this information we 
review the features of homogeneous and isotropic 3-d cosmologies, provide a number of new 
solutions and study gauge invariant perturbations around them. The first-order formalism is 
then used to provide a detailed study of the most general 3-d spacetimes containing perfect- 
fluid matter. Assuming the material content to be dust with comoving spatial 2-velocities, 
we find the general solution of the Einstein equations with non-zero (and zero) cosmological 
constant and generalise known solutions of Kriele and the 3-d counterparts of the Szekeres 
solutions. In the case of a non-comoving dust fluid we find the general solution in the case of 
one non-zero fluid velocity component. We consider the asymptotic behaviour of the families 
of 3-d cosmologies with rotation and shear and analyse their singular structure. We also 
provide the general solution for cosmologies with one spacelike Killing vector, find solutions 
for cosmologies containing scalar fields and identify all the PP-wave 2-1-1 spacetimes. 



1 Introduction 

General relativity in three spacetime dimensions is known to possess a number of special sim- 
plifying features: there are no gravitational waves, no black holes in the absence of a negative 
cosmological constant, the Weyl curvature is identically zero, and the weak-field limit of the 
theory does not correspond to Newtonian gravity in two space dimensions [l]-[9]. The theory 
is therefore considerably 'smaller' than general relativity spacetimes with four (or more) di- 
mensions, and the strong-energy condition that creates geodesic focussing does not depend on 
the density of the material sources. These simplifying features mean that considerable progress 
can be made in the search for the general cosmological solution of the three-dimensional Ein- 
stein equations. In an (A'^ -|- l)-dimensional spacetime the number of independently arbitrary 
A^-dimensional functions of the space coordinates that are needed to specify the Cauchy data 
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for the general cosmological problem on a spacclike hypersurface in vacuum is {N + 1){N — 2); 
in the presence of a general (non-comoving) perfect fluid it is N'^ — 1; and for a comoving perfect 
fluid it is - iV - 1 [2]. Thus, in the = 2 case, we see that the number reduces to zero for 
the vacuum solution (reflecting the absence of free gravitational fields in vacuum), reduces to 
one arbitrary spatial function in the comoving perfect-fiuid case ,and to three arbitrary spatial 
functions for a perfect fluid. 

In this paper we will set up the general cosmological problem in three-dimensional spacetimes 
and find the general solution of the field equations in the case of comoving pressure-free matter, 
with and without a cosmological constant, A. We go on to find solutions for the case of non- 
comoving dust and classify the singularities and asymptotic behaviours that arise in both cases 
with and without a cosmological constant. The relative tractability of the general cosmological 
problem in 2+1 dimensions allows us to go some way towards finding a general solution of the 
Einstein equations and we are able to isolate those features which prevent a full solution being 
found. In particular, we are able to find and classify the solutions for dust containing one of the 
(two possible) non-zero spatial 2-velocity components. 

There have been several past investigations of the structure of 2-1-1 dimensional general 
relativity and studies of the properties of particular solutions with high symmetry (see [1]- 
[8], and [10]- [14]). Important motivations for these studies were provided by the astrophysical 
interest in the possible observational signatures of cosmic strings and domain walls in the universe 
[15]- [18]. Higher-order curvature contributions were discussed in [2], together with the special 
features of the Newtonian-relativistic correspondence in general relativity and related theories, 
while the study of quantum gravity is reviewed in [19]. Cosmological solutions and singularities 
were discussed in [2] and [20, 21], while gravitational collapse of spherically symmetric dust 
clouds have been considered in [22]- [24]. 

The outline of this paper is as follows. In section 2 we define the 3-d Einstein equations and 
our notations. Section 3 introduces the 2+1 covariant formalism and the general kinematics of 
3-d spacetimes, identifying the special features that arise from the lower dimensions and from 
the vanishing of the Weyl curvature. These include the key role of the isotropic pressure as 
the sole contributor to the gravitational mass of the system and the fact that vorticity never 
increases with time. In section 4 we give a number of new cosmological solutions, review the 
characteristics of the homogeneous and isotropic models, including those that are singularity- 
free, and provide the generalisation of the Godel universe to 3 dimensions. We also consider linear 
perturbations around the 3-d analogues of the 'dust '-dominated FRW models and find them to 
be (neutrally) stable. In section 5 we employ Witten's first-order formalism [25] to formulate the 
equations governing the most general 3-d cosmological spacctimc metric containing perfect-fluid 
matter. Then, in section 6 we specialise the matter source to pressure-free dust with non-zero 
A and comoving 2- velocities and find the general solution of the field equations. These fall into 
three classes, one of which generalises the solution of Kriele [22] to A 7^ 0, while another is 
the generalisation of the Szckcrcs metric with nonzero A to 2+1 dimensions [26, 27]. Section 
7 considers the most general dust cosmologies with non-comoving velocities and finds various 
new classes of solutions. We study the asymptotic behaviour of these solutions and analyse in 
detail the structure of their spacelike and timelike singularities. Also, by means of a number of 
examples, we illustrate the wide range of possible behaviours in the presence of vorticity and 
shear. The same section also introduces a transformation that generates exact solutions with 
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nonzero cosmological constant from those with vanishing A. Finally, in section 8 we look at the 
case of a pure scalar field, provide the general solution of Einstein's equations with one spacelike 
Killing vector, and identify all the 2+1 PP-wave spacetimes. Our results are summarised and 
discussed in section 9. 

2 Einstein's equations 

It has long been known that general relativity, as a theory of a Riemannian spacetime, can be 
based on a small number of generally accepted postulates which are independent of the spacetime 
dimensions [28]- [30]. Assuming non-zero cosmological constant, these postulates demand that 
the field equations take the form 

Rab - \Rgab + ^gab = KTab , (1) 

where Rab is the Ricci tensor, with R = Ra"', Tab is the energy-momentum tensor of the matter 
generating the metric field gab, A is the cosmological constant and k is a dimensional coupling 
constant.^ When dealing with 3-dimensional spacetimes Qabg"^ = ^a^ = 3. In this geometrical 
environment R = —2kT + 6A, with T = Ta^, and Einstein's equations become (e.g. see [1, 2]) 

Rab = <Tab - Tgab) + Ikgab ■ (2) 

The spacetime geometry is determined by the Riemann curvature tensor Robed- In three 
dimensions the latter has six independent components, exactly as many as the associated Ricci 
tensor. This means that the spacetime geometry can be expressed solely in terms of the Ricci 
curvature, namely that [31, 32] 

Rabcd = gacRbd + gbdRac — gbcRad — gadRbc — \R{gacgbd — gadgbc) ■ (3) 

As a result, the 3-dimensional Weyl tensor vanishes identically and the gravitational field has no 
dynamical degrees of freedom. The spacetime curvature is completely determined by the local 
matter distribution and the theory is Machian. 

3 Covariant decomposition 
3.1 Observers 

In analogy to the standard 3+1 covariant approach to general relativity introduced by Ehlers 
[33] and elaborated by Ellis (e.g. see [34] for a recent review) , we introduce a family of timelike 
(fundamental) observers with worldlines tangent to the 3- velocity field Ua- The latter determines 
the time direction and is normalised so that Uqu" = — 1. The 2-dimensional space is defined by 
projecting orthogonal to Ua by means of the projection tensor 

Kb = gab + UaUb , (4) 

^Throughout this paper Latin indices run between and 2 and Greek take the values 1 and 2. Also, the three- 
dimensional coupling constant k is measured in units of mass~^ and therefore defines a natural mass unit [1]. 
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where hab = h(^ah)i habU^ = 0, habhf'c = ^ac and ha"" = 2. Using hab one also defines the covariant 
derivative operator of the 2-d space as = hJ'Vb- 

The irreducible kinematical variables, which describe the motion of the above defined ob- 
servers in a invariant way, are obtained by decomposing the covariant derivative of the 3- velocity. 
This splitting gives 

VfoUa = aab + ^ah + k^Kb - UaUb , (5) 

where aab = ^(b'^a) shear, ujab = Dj^u^j is the vorticity, O = D"i(a = V^Ua is the area 

expansion (or contraction) scalar and iia = u^Vb^a the 3-acceleration. Therefore, aabU^ = = 
(^abu'^ = Uau"" by construction. 

The tensor D(,Ua = hb'^ha'^VdUc = aab + '^afe + i&/'2)hab describes changes in the relative 
position of the worldlines of two neighbouring observers. When the latter follow the motion of 
a fluid, the effect of Q is to change the area of a given fluid element, without causing rotation 
or shape distortion. This scalar also defines the average scale factor, a, by 

- = ^©- (6) 

The shear monitors distortions in the element's shape that leave the area unaffected, while uJab 
describes changes in its orientation under constant area and shape. The symmetric and trace- 
free nature of aab ensures that it has only two independent components, while the antisymmetry 
of LUab guarantees that the vorticity tensor is determined by a single component. In other words, 
the shear and the vorticity correspond to a vector and a scalar respectively. The latter reflects 
the fact that the rotational axis has been reduced to a point. Defining eab = e[a6] as the 2- 
dimensional permutation tensor, with eabU^ = 0, the vorticity scalar is 

U = leabLO^^' , (7) 

with LUah = ^^ab- Note that = VabcU'^ by definition, where rjabc is the 3-d totally antisymmetric 
alternating tensor. The latter satisfies the condition r)abcV'^'^^ = ^^■^[a'^^b'^Scf : which ensures that 
eabe'"^ = 2hya''hf. 



3.2 Matter fields 

Suppose that the matter that sources the 3-dimensional metric field is a perfect fluid. Then, 
relative to an observer moving with 3- velocity Uq, the energy-momentum tensor of the material 
component takes the form 

Tab = PUaUb + phab , (8) 

where p is the energy density, p is pressure and its trace is T = 2p — p. Substituting the above 
into the Einstein field equations (2) the latter read 

Rab = 2{kp - A)uaUb + [k{ij. - p) + 2A]/l„6 , (9) 

with trace R = 2[k(/x — 2p) + 3A]. The above also provides the following auxiliary relations 

RabU^'u'^ = 2{kp - A) , ha'^hb'^Rcd = Hp -p) + 2N\hab and ha^RbcU^ = (10) 
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which win prove useful later. 

The twice-contracted Bianchi identities imply that V^Gab = and consequently the condition 
V^Tab = 0. The timelike and spacelike parts of the latter lead to the 3-d fluid conservation laws. 
These are 

P = -e{p+p), (11) 

for the energy density of the fluid, and 

{p + p)Ua = -^aP (12) 

for its momentum density. The above ensure that the conservation laws of a perfect fluid have 
the same functional form as their 4-dimensional counterparts (compare to Eqs. (37), (38) of [34]). 

The nature of the medium is determined by its equation of state. Here we will only consider 
barotropic fluids with p = wp, where w represents the barotropic index. When w = we 
are dealing with pressure-free dust, while isotropic radiation has p = p/2 and corresponds to 
w = 1/2 [1]. 

3.3 Spatial curvature 

The intrinsic curvature of the 2-dimensional space orthogonal to Ua is determined by the asso- 
ciated Riemann tensor. In analogy with its standard 3-d counterpart (see Eq. (77) in [35]), the 
latter is defined by 

T^abcd = ha'^hbhJh/Rqsfp - VacVbd + VadVbc , (13) 

where 

Vab = Dftlta = aab + tOab + ^©^afe , (14) 

is the relative position vector. Note that Vab characterises the extrinsic curvature (i.e. the second 
fundamental form) of the space. 

Starting from Eq. (13), assuming perfect-fluid matter and using expressions (3), (9), (10) 
and (14), the Riemann tensor of the 2-d (spatial) sections reads 

T^abcd = (^P - 3©^ + A) [hachd - hadhbc) - (o^ac + '^oc)(<7"&d + ^bd) + {(^ad + ^ad){(^bc + ^^bc) 
— \Q [((Tac + i^ac)hbd + hacicTbd + ^bd) — {o'ad + <^ad)hbc — had{<^bc + <^bc)] , (15) 

with IZabcd = T^[ab][cd]- agreement with standard 3-1-1 gravity, the isotropic pressure does not 
contribute to the curvature of the space orthogonal to Ua- Also, in the absence of anisotropy 
(i.e. when aab f'^id ujab vanish) the above reduces to 

T^abcd = {l^P - 3©^ + A) {hachd - hadhbc) ■ (16) 

Defining TZab = T^'^acb as our local 2-D Ricci tensor, we may contract expression (15) to 
obtain the following 3-dimensional analogue of the Gauss-Codacci formula (see Eq. (54) in [34]) 

nab = {f^-p-le^ + cT^ -u;^ + f^)hab, (17) 

which here holds for perfect-fluid matter. In deriving the above we used the results io^^^a^i^^ = 
and (Tc{a(^^b) = = ujc{a'^'^b)- The former holds because wi2(o"^i + 17^2) = (i.e. the single 
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independent component vanishes due to the trace- free nature of the shear). Similarly, the two 
independent components of crc{a'^'^b) ^Iso identically zero. Last, the result ojf.^^cij'^^^ = is 
guaranteed by the relation Uab = ^^ob and the properties of Eab (see section 3.1). The absence 
of a skew and also of a symmetric and trace-free part from Eq. (17) agrees with symmetries of 
the Riemann tensor in 2-d spaces (e.g. see [1]). 

Finally, the trace of (17) leads to the curvature scalar of the spatial sections, which may 
also be seen as the generalised Friedmann equation for three-dimensional spacetimes (compare 
to Eq. (55) of [34]) 

7^ = 7^«„ = 2 {np - + (7^ - + A) . (18) 

3.4 Kinematics 

The functional form of the Ricci identity is independent of dimension. Thus, when applied to 
the 3- velocity vector Ua, the Ricci identity reads 

VaVfetXc - VfeVaUc = Rabcdy!^ , (19) 

where Rabcd is the Riemann tensor of the 3-d spacetime (see expression (3)). Contracting the 
above along Ub, employing Eqs. (3), (10), and then taking the trace of the resulting expression 
we obtain the 3-d analogue of Raychaudhuri's equation 

e = -ie^ - 2kp - 2{a'^ - u;2) + B^iia + ifiia + 2A , (20) 

with 2(7^ = (TaftO""*' and 2a.-'^ = uj,jb^°'^- The above is the key equation of gravitational attraction, 
as it monitors the average separation between neighbouring particle worldlines. The most im- 
portant difference between (20) and its 4-d counterpart (see Eq. (29) in [34]) is that here only 
the fluid pressure contributes to the gravitational mass of the medium: the density p does not 
contribute. This unusual feature consists a major departure from standard gravity. One con- 
sequence is the existence of homogeneous and isotropic static 3-dimensional models with dust 
and zero cosmological constant (see solution (29) below). 

The symmetric and trace-free part of the contracted Ricci identity, together with the auxiliary 
relations (10) leads to the propagation formula for the shear in three dimensions: 

ha^h/acd = -Qo-ab + D((,n„) + in^Jib) ■ (21) 

Relative to the 4-dimensional case (see Eq. (30) in [34]), we notice the absence of the electric 
Weyl component from the right-hand side of this formula. This reflects the vanishing of the free 
gravitational field in 3-d gravity. Contracting this expression with Uab we obtain the propagation 
formula for the shear magnitude 

(a^)' = -29^2 + a'^'^B^^iia) + ^(^ti;.)^^"'' , (22) 

where 2(7^ = cTaba"'' by definition. In the absence of fluid accelerations the 2nd and 3rd terms 

on the right-hand side vanish and the equation integrates to give o"^ oc a~^. 

Similarly, the contracted antisymmetric component of (19) gives the 3-d counterpart of the 
vorticity propagation equation. The latter reads 

ha^hb^UJcd = -Q^ab + ^[bUa] ■ (23) 
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Comparing to Eq. (31) of [34] we note the absence of the ujc[a(^'^b] term. This is so because 

^c[a'^'^b] = in 3-d (see also above). Also, contracted with ujab^ recalling that 2lo'^ = u>ab^"'^ 
Eq. (23) gives the scalar vorticity propagation formula: 

(a;2)- =-2ea;2 + u;«fD[6U„]. (24) 

Again, in the absence of accelerations when the pressure vanishes this equation integrates to 
give cj^ oc as expected by the conservation of angular momentum. For a barotropic perfect 
fluid with p = p{p) and a constant barotropic index w = p/p, Eqs. (11), (12) combine with the 
commutation law Dj^D^]/ = —fuab to recast the vorticity propagation formula as 

K'ht,''u;cd = -{i-c!)@u;ab, (25) 

where Cg = dp/ dp = w is the square of the adiabatic sound speed. Therefore, the expansion 
decreases the vorticity of a 2-dimensional space unless the barotropic fluid has a stiff equation 
of state (i.e. when = 1 = w) and for general perfect fluids we have u>^ oc a'''^^^~'^\ Recall 
that, in standard general relativity, vorticity increases when ti; > 2/3 [36]. 

The simultaneous effects of shear and vorticity on a cosmology with negligible accelerations 
can be evaluated from these simple relations. For all equations of state we have cr^ oc but 
the centrifugal energy depends on the equation of state since oc a~^^^~'^\ Hence, the ratio 
of the distortion energy density to the centrifugal energy density is 

^ oc 

and the shear always dominates as a ^ when p > but the vorticity always dominates as 
a ^ DO. The presence of fluid acceleration can modify this behaviour; an example with p = 
will be given below. 



4 Cosmology in 2+1 dimensional spacetimes 

4.1 Homogeneous and isotropic spacetimes 

Spatial homogeneity and isotropy means that we can always choose a coordinate system such as 
the line element of the 3-dimensional spacetime takes the Friedmann-Robertson- Walker (FRW) 
form 

ds^ = -dt^ + a^hapdx'^dx'^ , (26) 

where the scale factor a = a{t) completely determines the time evolution of the model. The form 
also represents the metric of the 3-dimensional analogue of the familiar FRW cosmologies. The 
kinematics of (26) is monitored via one propagation and one constraint equation (see expressions 
(20) and (18)). Assuming a barotropic fluid with with constant barotropic index, setting A = 
and recalling that 9 = 2d/a, these formulae recast into 

a „ , f k 

- = —ZKwp and - = Kp , (27) 

a \a J a'^ 
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respectively {k = 0, ±1 is the curvature index of the spatial sections). The matter component 
obeys the conservation law (11), which accepts the solution 

P = m{i) . (28) 

with oq constant. 

For dust, w = 0, and one immediately finds (see Eq. (27a)) that a oc i for all signs of k. 
More specifically, expression (28) gives p oc a~^, which substituted into Eq. (27b) leads to 

a = [kpoUq — fc)^^^ t + ao , (29) 

where the product poo^ is proportional to the total mass of the model. Therefore, the scale factor 
of a dust-dominated, 3-dimensional, FRW universe evolves linearly with time, irrespective of its 
spatial curvature and no collapse to a final singularity occurs. Note that when kpoQq = k we 
obtain a static solution with a = ao [1]. Unlike its 4-d counterparts, this static universe has zero 
cosmological constant (see also section 5 above). 

When the matter content is in the form of black-body radiation, the barotropic index is 
w = 1/2. Then, expression (28) gives p oc a^^ and the system (27) has the solution 



- = J f''l~^" dt + -o. (30) 

For k = the above reduces to a oc i^^^, which coincides with the scale-factor evolution in 
standard 3-1-1 dust-dominated FRW universes. 

One can study perturbations around the above given homogeneous and isotropic solutions 
by introducing the covariantly defined variables T>a = {a/ p)T)ap and Za = aDa©. The former 
describes variations in the matter density and the latter in the area expansion as measured by 
a pair of neighbouring observers [37]. Also, both variables vanish in the spatially homogeneous 
background and therefore satisfy the gauge-invariance criterion. In the case of a pressureless 
fluid (i.e. for u; = = Cg) the linear evolution of inhomogeneities is monitored by 

V = -Z and Z = -@Z, (31) 

on all scales. Then, using solution (29) and setting Pq = ("^o + aQ)ZQ/a initially, we find that 
the density gradient decays as 

P = Pof^^), (32) 

V at + ao / 



with a = ■\/ kpqo^ — k and k = 0, ±1 (see (29)). For general initial conditions, on the other hand, 
it is straightforward to show that P — > Pq — {ato + ao)Zo/a at late times. Recalling that for zero 
pressure shear and vorticity perturbations decay in time (see Eqs. (21) and (23)), we conclude 
that in the absence of pressure the 3-d analogues of the FRW cosmologies are either stable or 
neutrally stable. This behaviour is very different from that of conventional FRW models, all of 
which are unstable to density perturbations, and reflects the fact that in three dimensions the 
gravitational mass of a pressure- free medium vanishes (see Eq. (20)) and spherical regions of all 
curvatures asymptote to a ^ t as t ^ oo. The immediate consequence is the absence of linear 
Jeans-type instabilities in these models. 



8 



4.2 Homogeneous and anisotropic spacetimes 

The simplest 3-d line element describing a homogeneous and anisotropic spacetime has the 
following Bianchi I-type form 

ds^ = -dt^ + A'^dx'^ + B^dy^ , (33) 

where A = A{t) and B = B{t) are the two individual scale factors. When A = 0, the spatial 
flatness of the above metric means that the 3-d analogue of the Bianchi I universe is covariantly 
described by the following set of propagation equations 

p = -{l + w)Qp, e = -^e'^ -2Kwp-2a'^ , a = -ea, lo = (34) 

supplemented by the constraint 

le^=Kp + a\ (35) 

The above are obtained from Eqs. (11), (20), (22) and (18), respectively, after dropping their 
inhomogeneous terms and assuming spatial flatness together with zero vorticity. Setting O = 
2a /a, where a represents the geometric-mean scale factor, expressions (34a) and (34c) lead to 
the evolution laws 

/ao\2(i+«') /ao\2 
P = Po y—j and cr = (To y— j , (36) 

for the matter energy density and the shear respectively. Substituting these results to Eq. (35) 
and assuming a spacetime filled with pressure-free dust we arrive at 



^Kpof^ + 2aot . (37) 

As expected, for (Jq = the above result reduces to its isotropic counterpart (see Eq. (29)). 
The same also happens at late times, when the shear contribution to solution (37) becomes 
negligible. 

Note that one can use the evolution law of the average scale factor to obtain those of the 
individual ones. Assuming dust, recalling that 

^ + | = e and ^-| = 2c.> (38) 

and using results (36c) and (37) we arrive at 

t\ J „ „ / npot + 2(70 

'0 r 



A = A„^-^ and B^B„y.J^^^y (39) 

In the vacuum case (po = 0) the metric is just a coordinate transformation of fiat spacetime 
rather than an anisotropic cosmology. 
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4.3 Rotating spacetimes 

Consider a rotating three-dimensional spacetime with flat 2-d sections filled with pressureless 
matter. Setting A = and assuming spatial homogeneity, the time evolution of the model (at 
least locally) is monitored by the following system of four propagation equations 

p = -ep, e = -ie^-2c72 + 2u;^ a = -@a and uj = -Qu , (40) 

constrained by 

^e^ = Kp + (7^-a;^ (41) 

Proceeding as before, we may use the area expansion scalar (9) to define an average scale factor 
(a) so that 9 = 2d/a. Then, expressions (40a), (40b) and (40c) translate into 

P = Po (^^) , (T = ao and cv = uq (^^^ (42) 

respectively. Substituted into constraint (41) the above results give 

a = aQ\j Kpof^ + 2t\J al-ul- (43) 

Accordingly, despite the presence of nonzero shear and vorticity, the average scale factor evolves 
as its homogeneous and isotropic counterpart if uq = uq. The situation in 3-d is unusual in 
that the shear and vorticity scale as the same powers of the scale factor and are both equally 
important at all times. In the case of dust, the matter density also scales as a^^ (see (42)). 

One can also obtain a 3-d Godel-type universe. The Godel spacetime is a homogeneous space- 
time and a rotating solution of Einstein's equations which permits closed timelike curves [38, 39]. 
Covariantly, Godel's world is described by [40, 41] 

Q = {) = Ua = Gab and cjab 7^ . (44) 

Thus, with the exception of the vorticity, all the kinematical variables vanish identically. Note 
that the overall homogeneity of the Godel spacetime guarantees that the vorticity is a covariantly 
constant quantity and that all the propagation equations reduce to constraints. Applying (44) 
to three dimensions we arrive at the system 

p = 0, d; = 0, Kp-J^ - K = and 7^ = 2(fcp - u? + K) . (45) 

which describes a 3-d Godel-type universe. Note that for pressure-free matter the cosmologi- 
cal constant is necessarily negative (i.e. A = — - see Eq. (45c)). Also, unlike its standard 
3-|-l counterpart, the 3-dimensional Godel-type spacetime can have non- vanishing spatial Ricci 
curvature. This is guaranteed by (45d). The empty Godel-type model with TZ = — 4c<j^ < is 
equivalent to anti-de Sitter (AdS) space. Rooman and Spindel, [42], considered a one-parameter 
subset of these non-flat Godel-type solutions and showed that they can be seen as arising from 
a directional squashing of the light cones of AdS, which breaks the so(2, 2) isomctry of the AdS 
Killing vectors into so (2, 1) x so (2). It was shown in [42] that all the non-empty Godel-type 
solutions considered contained closed-timelike curves but they vanish in the AdS limit. 
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4.4 SingulEirities 

The 3-d analogues of the standard singularity theorems are relatively straightforward to deduce. 
The study of the Riemann and the Ricci curvature shows that there arc no Weyl curvature 
singularities and the analogue of the strong-energy condition (i.e. Rabu""u^ > 0) reduces to the 
inequality p > for a perfect fluid and to the positivity of the sum of the principal pressures 
if they are anisotropic. Also, the form of the Raychaudhuri equation (see (20)) guarantees that 
for geodesically moving observers 

Q + ^Q^ = -2kp-2{(j'^ -uj'^) + 2k. (46) 

Therefore, for vanishing cosmological constant and in the absence of rotation, an initially con- 
verging family of timelike worldlines will focus (i.e. O — oo). For non-zero vorticity, however, 
this may not be necessarily the case. For example, applied to a spatially homogeneous and 
rotating spacetime filled with an isotropic radiation fluid, the above gives 

G + iG^ = -.po(^)'-2(ag-.o^)(^)'. (47) 

Therefore, as a — > 0, the kinematical term dominates the matter term in the right-hand side 
of the Raychaudhuri equation. In this case, whether caustics will form or not depends entirely 
on the balance between shear and rotation. When ujq > (Tq, in particular, vorticity can stop an 
initially converging family of worldlines from focusing. 

There are also simple exact solutions that describe 'bouncing' 3-d cosmologies with p < 0. 
Consider, for example, a perfect fluid with p = —p/2. Then, assuming spatial homogeneity and 
isotropy, Eq. (11) gives p oc and the Friedmann equation reduces to 

M k 

= 2 ' (48) 

where M > is constant and k is the spatial curvature index. The bouncing solution is 

a{t) = a^in + lM{t- t^inf , (49) 

with flmin = k/M. Accordingly, a non-singular minimum requires positive curvature. Setting 
k = +1 and tmm = the required scale- factor evolution is represented by the parabola 

a{t) = M-^ + iMf . (50) 

5 First- Order Formalism 

In 2-1-1 dimensions it has been shown by Witten [25] that the Palatini action for general relativity 
is equivalent to: 

Sg= [ r'^ei^F,,! (51) 
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where 77"^'^ is the metric-independent alternating symbol in 2 + 1 dimensions. Recall, our 
signature is ( — h +) and wc choose fj^^"^ = — 1 =^ 77012 = +1; is a dreibein, and in some sense 
should be viewed as the square root of the metric 



gob = eie^ViJ, 



where 77/j = diag{—l + 1 + 1). We set 



where ^A^ is an 5*0(2,1) connection, and ejjx is the standard alternating symbol. Here 
o, 6 = 0, 1, 2 are spacetime indices, whilst I, J,K = 0, 1,2 arc internal 5*0(2, 1) indices. We raise 
and lower internal indices with 77/ j and r]^'^ , so e^"^^ = —ejjK, and we choose eQ^2 = +1- The 
action (51) reduces to the standard second-order gravitational action when is compatible 
with i.e. 

T^laeb]! •■= 5[ae6]/ + eiJK 'Af„ef = 0. (52) 

We could also find this equation by varying the above action w.r.t. ^A^^. There is plenty of 
gauge freedom in this model, so we firstly simplify matters by choosing to write the metric in 
synchronous gauge as 

ds^ = gabdx'^dx^ = -dt^ + ha^dx^'dx^, (53) 

where a, P = 1,2. In synchronous gauge: Cq = 1, Cq = 0, e[J^ = 0, where i = 1,2 are 2-d internal 
indices. Throughout this section, we will use a, b, c, d to designate 3-d spacetime indices, and 
/, J, K for 3-d internal indices. Lower case Greek letters will be used for 2-d spacetime indices, 
and k for 2-d internal ones. We define a projection operator onto the surfaces t = const by 

Kb = gab + naUb (54) 

where Ua = (1,0,0) and hob = h^p is the induced metric on the spacelike surfaces defined by 
t = const. We make some further definitions: 

Ei = h%, (55) 
Al = h^jA^^, (56) 

Ba = h^jApo^ (57) 
= hlh^F^si, (58) 
{n-AY = n^'^Al (59) 

We can also use the gauge freedom in the definition of ^A^^ to set {n ■ = 0, and do so. We 

are still left with the freedom to make 50(2) rotations on the internal indices, z = 1, 2, with the 
angle of rotation an arbitrary function of the spatial coordinates, {x,y). 
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5.1 Field Equations 

In 2+1 spacetime dimensions: 



Rac9bd — Rad9bc + Rbdda 



Rbc9ad + ( "2 ) (dadgbc — gac9bd) , 



Rab = l^(Tab — Tgab), 

R = -2kT. 

So, therefore: 

Rabcd = 1^ [Tac9bd - Tad9bc + ^bd^ac - ^bc^ad + T {gad9bc - 9ac9bd)] : 

where the curvature ^Fabi is related to the Rabcd by 

^Fabl = Rabcde'^'^e'^^eijK- 

We define: T^"^ = e""^ e^^Tab, so the Einstein equations can then be written as: 



(60) 

(61) 
(62) 



^[066]/ = 0, 

^Fabl = K 



(63) 
(64) 



Note that these equations remain well-defined even if detie) = \/—g = 0, so in this form the 
field equations actually describe a theory that is more general than general relativity. When 
is invertible, however, these field equations are fully equivalent to 3-d general relativity. We 
are concerned with the dynamics of perfect-fluid spacetimes and for the purposes of this section 
will consider only dust models, p = 0, so T^^ = npu^w^ — Arj^-^ , where uku^ = —1- We write 



u 



= Uk and Ua = l + UkU\ k = 1, 2. 



5.2 2+1 decomposition of field equations 

We now use n" and to decompose the field equations. Prom equation (63) we find that 
n- = and: 



E ■ — f - A^ 



r^daEpi + eijB^E^^r^ = 0, 
V^^eijAiEj, = 0, 

and by projecting equation (64) we arrive at: 

Ba = KpuQWeijEl, 



r^{2d^Bp + e,,A'^£,) 



Kp 



(65) 
(66) 
(67) 

(68) 
(69) 
(70) 

lElE'^r^UkU'eij + E^E^^r^ei^ + Ae^.E^E'^r^ . (71) 



np UkUkeijE^p - UiE^p€ijU^ -KeijE^p, 
-Kpu^U^ejkE^Epir^ 
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Equations (65), (68) and (69) provide 10 evolution equations for the ten degrees of freedom 
contained in the E^^, and A^. To solve these we must specify these quantities, as well as 
p and the C/*, on some initial spacelike hypersurface: a total of 13 initial functions. The other 
equations amount to 6 consistency conditions, which restricts the number of free functions to 7. 
We may still make 3 coordinate transforms, and 1 internal SO{2) rotation, or gauge transform. 
This brings us down to 3 functions that may be freely specified on the initial surface; we take 
these to be p, and U"^. It is easy to check the the consistency conditions are preserved by 
the evolution equations. 

6 General Cosmological Solutions with Comoving Dust 

The equations are particularly simple when the velocity can be chosen to be co-moving (which, 
for 2+1 dimensional dust, is equivalent to it being irrotational) . We fix our coordinate system 
by choosing both hag and hajs to be diagonal on our initial surface; for a proof that this may 
always be done see [22]. It was also proved by Kriele that if, initially, hap 7^ X{x,y)ha/3 then 
this coordinate choice is unique and the only remaining freedom is t ^ t + to, x ^ ^i^) and 
y — > Y{y) although we also retain the freedom to interchange x and y. In this case, and only 
this case, we may use the coordinate and gauge freedoms to move to a frame where the fluid is 
comoving and i?* and El^ are diagonal at some initial instance, and then the evolution equations 
ensure that they remain diagonal at all times. Such a coordinate transform leads to eqn. (67) 
being automatically satisfied. Thus, in making this coordinate choice, we reduce the number of 
consistency equations by 1 and so the number of free-functions that can be specified on some 
initial surface increases by 1. We have already set = U'^ = so we are therefore left with 
2 free functions. Solutions of this system represent all 2 -|- 1 dimensional dust spacetimes with 
vanishing vorticity. 

If hapito) = X{x,y)haf3{to) for some X{x,y), then we may further transform to a conformal 
gauge so that 

hag = (1 + X{x,y){t - to))e''^(^'^)(5a/j, 

where dap = diag(l, 1); such a coordinate choice is unique up to a; — > Ax + B, y ^ Ay + C, 
t ^ t + to- We keep the freedom to interchange x and y. We shall see that such spacetimes are 
specified by only one free function of x and y. 

We use our remaining coordinate freedom to set h^g and hag to be diagonal initially, and 
fix the gauge of our internal index by requiring E'j^ and eijAa = E^i to be diagonal initially. For 
comoving systems J7* = 0, and so by combining equations (65) and (69) we see that 

Thus, if E'q and CijAa = Eai are diagonalised initially, they will remain diagonal for all time. 
We have now fixed our system and our internal gauge freedom. With these choices equation 
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(67) is automatically satisfied. Solving for E'^ gives us 

eI = C(x,2/)shA(t) + D(x,y)chA(i), (72) 

eI = W{x,y)shA{t) + V{x,y)chA{t), (73) 

4 = Cix,y)diAit) + ADix,y)shAit), (74) 

4 = -W{x,y)chA{t)-AVix,y)shAit). (75) 



with C{x,y), D{x,y), W{x,y) and V{x,y) to be determined by the remaining equations. We 
have defined shA(t) = A~^/^ sinh(A^/^t), and chA(t) = cosh(A^/^t). All other components of 
and El^ equal to zero. For comoving systems we also have 

Ba = 0. 

This equation is automatically satisfied whenever the two consistency conditions, (66) and (70), 
hold. In our choice of coordinates and with our gauge-fixing these equations read: 



dyC{x,y) = BiW{x,y), (76) 

dyD{x,y) = BiV{x,y), (77) 

d^W{x,y) = -B2C{x,y), (78) 

d^V{x,y) = -B2D{x,y). (79) 



The 2-surface, t = const, is flat whenever fj"^daBp = 0. The solutions of these equations divide 
into 3 distinct classes. 



6.1 Class 1: = and / or = 

If one or other of the Ba vanishes we can, without loss of generality, by interchange of x and 
y, choose Bi = and leave B2 to be freely specified (we are free to choose ^2 = if we wish). 
Except when B2 = B2{x), these solutions will generically require the 2-surface given by t = const 
to be non-flat. The solution of equations (76)-(79) is simple in this case. We find C = C{x), 
D = D{x) and: 

W{x,y) = £d^CiO^^^ + f{y), (80) 

V{x,y) = 1^^0(0^^^ + g{y), (81) 
_ dL{x,y) 

^^^^ 

We determine the energy density of the dust, from eqn. (71), to be 

_ -dlL{x,y) + (l^(a;,|/)chA(t) + Ay(x,j/)chA(t))(C(x)chA(t) + AD(x)shA(t) 
(C(x)shA(i)+i^(x)chA(t))(iy(x,y)shA(i) + y(x,2/)chA(0) 

and the metric is 

ds2 = -dt^ + (C(x)shA(t) + D{x)c}iK{t)fdx'^ + {W{x, y)shA(t) + V{x, y)c\iA{t)fdy'^ . 

This is equivalent to the dVo/dy = solution found by Kriele, [22]. As stated above, the limit 
i?2 — * is well-defined in this case. 
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6.2 Class 2: Si 0, 7^ and C /W ^ D /V 

It will be seen that Class 1 solutions do not emerge as the 5i = limit of the solution for Bi and 

B2 non-zero. The solution of the latter case is more complicated. The condition C /W 7^ D/V 
ensures that we do not have hap oc hap-, and so our coordinate choice is unique (up to rescalings 
of X and y). 

As discussed above, solutions of this class are defined by 2 free functions of {x,y)\ to aid 
finding the solution we choose these to be C (x, y) and D(x, y). By combining equations (76)-(79), 
we can solve for Bi: We solve this equation for Bi: 

where S{y) and xq are arbitrary. Equations (76) and (77) can now be used to specify W and V: 

W = e-%C/S{y), (83) 
V = e-%D/S{y). (84) 

Finally, we check that these forms of C, D, W, and V define a unique B2 via equations (78) and 
(79). We find that they do, and that B2 is given by 

e"^ (dyCd^dyD-dyDd^dyC 

-D2 = 



S{y) V CdyD - DdyC 

For D 7^ 0, we may write these in a slightly more familiar form by defining D = e^^^'^^ and 
C = F(x,y)e^(^'f). We then have 



= lnF.y + j/ + a, (85) 
a{x,y) : = Td^^^. (86) 



Thus, Bi = S{y)Fye''+'' and 



and 



^ e-'^{Fen,ye-" 
F,yS{y) ■ 

FySiyY 



_ e—-(y,yf f Fy 



The metric is therefore given by: 

S^y)Fl{x,y) 



d.^ = -dt^ + R\x, y, t)e-(-^)dx^ + ^"''""'ll^^l^^f/^^^'^'^^-^^d/ 



where R{x, y, t) = ch\t + F{x, y)sh.\t, and the functions F{x, y) and i/(,t, y) arc arbitrary. Gen- 
eralised Szekeres-like solutions emerge in the limit F^x = 0. The requirement C/W 7^ D/V 
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translates to e'^ ^ H{x){F{x,y) — 1) for some H{x). If this requirement does not hold then, 
although our solution is still valid, the coordinate system is not uniquely specified and we may 
transform to a conformal frame in which the solution is simplified; we deem such solutions to 
be of Class 3 and deal with them in the following subsection. The energy density for the Class 
2 solutions is now given by a rather complicated expression: 

R{R,y + Rl',y) ' 

where 

E{x, y, A) = e"^'' ['^,xyi^,x - v,xxy - Fx/Fl{iy^xi^fyF^y + F^^u^y - Su^^yi^^yF^y) (88) 

~ E,xy/ Fy{y^yV^xF,y + ^^,yE,x + '^F,xy^,y ~ ^^,xyF,y) 

+ F,XX^l/F,y + F^xxyU^y/F^y] + ^C'^'' [K C^'') ^ 

with K{x,y,t) = B? — h.R^ — e^^{SF^y)'^. Generalised Szekeres solutions [26, 27] emerge in the 
limit Fa; = and the 2+1 Szekeres solutions [43] themselves emerge when = A{y)x^ + 
2B{y)x + C{y). 



6.3 Class 3: Bij^0,B2^0 and C/W = D/V 

In this case hai3 oc hai3 and so we may, as mentioned above, transform to a conformal gauge 
where both the metric and its time derivative are proportional to 5ab- In this case, without loss 
of generality, D = V = e^^^'^^ and C = W = ne'^^^'^^ Bi = dy(p and B2 = —dx4>', <P is arbitrary 
and /X is a constant. The metric then becomes 

ds2 = _di2 + e2*(MshA(i) + chA(t))^ (dx^ + dy^) 

and the energy density is: 

_ -e-2^V20(x, y) + (/xchA(t) + AshA(t))2 
(MshA(i) + chA(i))' 

where = 5^ + 5^. The FRW solution emerges from the case where e '^^V'^(t){x, y) = const. 
This is also the only homogeneous limit of this class of solutions. 



7 Cosmological Solutions with Non-Comoving Dust 

In this section we will extend our discussion to include the remaining degrees of freedom in the 
cosmological evolution of dust. We consider non-comoving fluid motions and seek a new class of 

non-comoving solutions where only one of the spatial velocity components is non-zero; without 
loss of generality, we will take this to be the a;-component. Wc will initially work with A = 0, 
however we shall present a simple transformation that allows us to map A = solutions into 
A 7^ ones. 
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First, we use some of the remaining gauge and coordinate freedom to set E2 = 0. With this 
choice we have Ux = sinhS ^ 0, U2 = and ug = cosh0. This is equivalent to demanding that 
we choose our local frame field so that the j/-component of the fluid velocity vanishes. As before, 
we have that Eai = CijAa, and the relation ij^f^eijAl^Ej^ = tells us that = Ef/Ej, 

which implies E^ = A{x,y)E2. Hence, we find that El = C{x,y)t + D{x,y). The remaining 
(independent) equations to be solved now read: 

e| = -Kpsinh'^eEl, (89) 
B2 = Kp cosh 9 sinh 9 eI, (90) 
Si = A{x,y)Byix,y,t) + fix,y), (91) 
dy{Cix,y)t + D{x,y)) = f{x,y)El (92) 
(dx - A^y - Ady) eI = -B2iC{x,y)t + Dix,y)), (93) 
KpcoshH{C{x,y)t + D{x,y))El = {d^ - A^y - Ady) By - dyf{x,y) (94) 

+C{x,y)El. 

By taking two time derivatives of equation (92), we can see that we must have f{x,y) = 0, 
which in turn implies C = C{x) and D = D[x). This is therefore a generalisation of the Class 
1 comoving solutions found in the last section to case of non-zero vorticity. In the system of 
equations above, we have written down two evolution equations and two consistency equations, 
and by combining the evolution equation for B^ with the time derivative of the first of the 
consistency equations we arrive at a third condition. On any initial surface, we must specify 6 
functions of x and y: A{x,y), p, 9, £J| and and B2, and 2 functions of x: C{x) and D{x). 
Given the assumed form of the metric, the only coordinate freedom we have left is y ^ Y{x,y)^ 
X X{x) and t ^ t + Iq. At some initial instant wc can always use the residual coordinate 
freedom to set A{x, y) = 0. Without loss of generality, therefore, we set A{x, y) = 0. This leaves 
5 free functions and 3 consistency equations. In total, we are left with 2 free functions that can 
be freely specified on the initial surface, one less than required for the general solution of the 
Einstein equations in accord with our setting one of the 2-velocity components equal to zero. 

We shall define new variables: X = Hcpcosh^ 9, Y = tanh^, B2 = L(x,y,) and E2 = 
M[x,y,t). With these definitions the equations become: 

M = -XY'^M, (95) 

L = XYM, (96) 

dxM = -L(C{x)t + D{x)), (97) 

X{C{x)t + D{x))M = dxL + C{x)M. (98) 

When C 7^ we can, at least locally, set (7=1, without loss of generality, by using the freedom 
to redefine the x coordinate. In this case, we can combine the above equations (95)-(98 )into a 
single second-order, non-Hnear, PDE for M{x,y,t): 

dt (ff^) )\t + Df= Ut + D)dx (^) -{t + D)^,M^ dfM. (99) 
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We note that via the coordinate transformation {t,x) — ^ {t',x'), defined below, we can, without 
loss of generality, set D{x) = and preserve all the assummed properties of the metric and 
matter content. 

t'coshx' = «coshx + /"ci«D(0smh5, 

/X 
d^D{^) cosh^. 

Hereafter we fix our coordinate system by taking D = 0. 



7.1 General solution for dust with one non-comoving velocity 

By the coordinate transform defined above we set D = and solve the resulting system of 
equations. The only subcase not explicitely covered by this solution will be that when C = (in 
which case we can take D = 1 w.l.o.g.). We shall see later through the class of solutions with 
C 7^ is, up to a coordinate transform, equivalent to the class of solutions with C = 0, D = 1. 
We define a new variable Z = XM. The system of equations now reads: 

M = -ZY^, 
L = ZY, 
d^M = -Lt, 
tZ = dxL + M. 

By combining the first three equations we find: Y = u^t/u^x where t = \nt and u = M,^ — M. 
We then take the r-derivative of the last equation to arrive at: 

drn-Yd^n = d:,Y -Y'^, 

where we have defined = Zt. We now make a coordinate transform: {x,y) {u,X), where 
u = M^T — M and X = x. With respect to these new coordinates, the above equation becomes: 

In terms of these new coordinates 1/Y = — tx- We can rewrite M = —ZY"^ as u^-j-e"'^ = —ZtY^. 
Inserting this into the above equation leads to a simple equation for Y: 

Yx = Y^- 1, 

which has solution Y = tanh(a; — d{u,y)) with d{u,y) a function of integration. Solving the 
above equations we then find 

Zt 

— = — = C{u, y) cosh(x — d{u, y)), 

with C{u,y) a function of integration. Prom u^r^-"'^ = —ZtY^ we have 

e~'^ = ~C{u, y) sinh(x — d{u^ y)). 
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It will be more straightforward to define F'{u,y) := F„ = — C{u,y) cosh. d{u,y) and G'{u,y) := 
G,u = C{u,y) smb. d{u,y). It is clear that F'^ — G'^ > 0. With respect to these definitions the 
above equation for r becomes: 

e"^ = ^ = F' sinh(ar) + G' cosh(x). 

From this, we find expressions for u^r and u^x- 

F' sinh(x) + G' cosh(.T) 



F" sinh(x) + G" cosh(a:) 
F' cosh(x) + G' sinh(x) 



F" sinh(x) + G" cosh(x) ' 

Finally, we find M := tP by solving dT{e~'^ M) = u to obtain 

F G 
P = -u^{—)' sinhfx) - u^{—)' cosh(x) + H{x, y) 
u u 

where the function H{x, y) is found, by insertion of M into tZ = dxL + M, to satisfy H^^x = H. 
We can therefore absorb H{x,y) into the definition of F{u,y) and G{u,y), and without loss of 
generality set H = 0. Thus, we have the final form of the general non-comoving dust solution 
with Uy = : 

M = -tn^ ^^^j sinha;+ l^^j coshxj , (100) 

F'^ - G'^ 

'^^ ~ {F'^{u/Fy sinhx + G'^iu/G)' coshx) (F" sinhx + G" coshx) ' ^^^^^ 
y _ F' sinhx + cosh a: ^^^^ 
F' cosh X + G' sinh x ' 



^ sign(r) 



^ |F-coshx + G-sinhx| 

It is evident from the form of the density, p, that there is a curvature singularity at t = 0, and 
that to ensure we always have positive energy densities we must have 

{F'^{u/F)' sm\ix + G'^{u/G)' cosh x) (F"sinhx + G" coshx) > 0, 

with cTirvaturc singularities appearing in the case of equality. We could, it should be noted, 
rewrite this as the requirement that 



(p2y 

^—L = {u^{F/uy sinhx + u^{G/uy coshxf /u < 0, 
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with singularities forming in the case of equaHty. In the next subsection we shall consider the 
form and nature of these singularities in more detail. Using the above results we can find the 
expansion scalar, O, for these spacetimes: 

A{u,y) sinhx — B(u,y) cosh X 
tP{u,x,y)P' 

where 

A{u, y) = {2FG'F" - GF'F" - FF'G" + F'^G"u - G'F'F"u) , 
B{u, y) = {2GF'G" - FG'G" - GG'F" + G'^F"u - G'F'G"u) . 

We also find that the vorticity of these spacetime is: 

2 _ - G'F'F[yf - - F'G'Qyf 

7.2 Classification of Singuleirities 

In addition to the standard cosmological singularity we have singularities in the dust density 
(p — ^ oo) whenever: 

= 0. 

We can divide these (-P^)' type singularities into two distinct classes: class A singularities arc 
where P / 0, P' = 0, and class B are where P = 0. Class B singularities can be thought 
of as shell-focusing singularities, in analogy to those in the Szekeres spacetimes. For class A 
singularities, the volume element of the metric remains non-zero, whereas for class B singularities 
it vanishes. From the definition of u, we have that u = iff u^r = which will not be the case 
at any finite time. Thus, at finite times, either « > or ti < 0. 

Class A and B singularities are generically naked. This can be seen by consider geodesies 
that move along y = const paths. The metric along dy = is: 

ds^ = -dt^ + t^dx^ = -t'^dudv, 

where u = Int + x and v = Int — x. From any point, {to, xq, yo}, we can therefore move along 
an outward-moving null geodesic, defined by y = 0, u = const, that will reach null infinity. It 
follows that there exist no black-hole horizons in this spacetime. We now consider the strength 
of a singularity at the point {to,xo,yo} = {uo,xo,yo}. We consider the quantity 

* := RabK'^K^ 

where = dx"'{k)/dk, and k G (0, 1]. The singularity lies at A; = 0. Since the Weyl tensor 
vanishes in 2-1-1 dimensions, by propositions 1-4 of Clarke and Krolak, [44], a necessary and 
sufficient condition for the singularity to be strong in the sense of Krolak is that the following 
integral does not converge as — > 0: 

J{k) = J dk'^{k'). 
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If limfc_»o J does not exist then the limiting focusing condition (LFC) is said to apply. If J{k) 
does converge as A; ^ then the singularity is weak in the sense of Krolak, and also in the 
sense of Tipler [45]. A necessary and sufficient condition for the singularity to be strong in the 
sense of Tipler is that the strong LFC apply, i.e. that J{k) not be integrable in (0, 1]. If J{k) is 
integrable in (0, 1] then the singularity is Tipler weak. 

The equations describing null geodesies in this background are: 

= p K^ = f.^ (105) 
V = 2' + ^, (106) 



M2 

dV uP 



dk t^P^ ' 

dQ _ {F' cosh x + G' sinh x)f 

dfc ~ ' 
dl {F^y sinh x + G^y cosh x)f' 

dk ^ Wp^ ■ 



(107) 
(108) 
(109) 



We can see that simple solutions can be found if we take 1 = = . In these cases V = r^Q 
const. Therefore, along these null geodesies we have 



t = ytl + 2Pk, x = xo±^ln{l + 2Pk/tl). 

Outward-moving null geodesies take the + sign, while inward moving ones correspond to the — 
sign. We shall refer to these geodesies as 'radial' null geodesies (RNG). We define the quantity: 

A = (F'coshx + G'sinhx)(l±y). 

For finite t > 0, A is both finite and non-zero. Along RNGs we have: 

du PXu 

dk = ^-PF- 

^ = T^IA-F.). (Ill) 



dP' P 



- (XuP" - P'P[,) . (112) 



dk 

Proposition: The LFC does not apply to RNGs terminating on Class A singularities. 
Proof: Consider the quantity ^ for class A singularities: 

^ - ^'(-^)^' 



Using 1/P' = :3j:t^/{P\u){du/dk) we have: 

\P du 



* = ^ 



Pt2 dk' 
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The limit limk-^o XV / Pt^ exists for class A singularities and du/dk is integrable on (0, oo), and 
so ^ is integrable on the same region. Therefore by propositions 4 and 6 of [44] the LFC docs 
not apply to RNGs for class A singularities. Class A singularities are therefore gravitationally 
weak □. 



Proposition: The LFC applies along RNGs terminating on Class B singularities provided 
limfc^o l/-f" / 0; but the strong LFC does not apply. 
Proof: It is a sufficient condition for the LFC to hold that: 

hm k"^ = lim — 7^ 0. 

fe_0 fe-»o PP't"^ ^ 

We note that, for t G (0, 00): 

XqUq - Pa;|o = F{uQ,yQ) cosh, xq + G(«o,yo)sinha;o ±ito/io / 0, 

and is finite. By I'Hopital's rule and equation (111), the limit \\m^^^Qk/ P must therefore exist 
and be non-zero. Thus limfc^o = iff lim/j^o 1/-P' = or equivalently iff limfc_>o du/dk = 0. 
The LFC applies along all RNGs where \mik^Qdu/dk ^ 0. By propositions 1 and 2 of [44] the 
strong LFC does not apply if the integral 



J{k) ■= j\k' dk"^{k") 



converges as A: ^ 0. Using integration by parts we see that this is equivalent to the condition 
that /-^ dk'k''^{k') converges. Using 1/P' = ^t^ /{VXu) {du/dk) we have 

k XV du 

* = ic 

^P t2 dA;' 

As shown above, the limit limfe^o k/P exists and du/dk is integrable on (0, 1], therefore k^{k) 
is integrable on (0, 1] and the J{k) converges, and the strong LFC does not apply □. 

7.3 Asymptotic Behaviour 

We consider next the asymptotics of our class of non-comoving dust cosmologies and find the 
criteria for them to become homogeneous at late times. Wc note that as t ^ 00, we have 
F' sinhx + G' cosh re 0"*". We assume that as t — > 00, the quantity F" sinhx + G" coshx does 
not vanish. We also assume that the limit lim^^oo^* = 1*0 exists and that all functions of u 
have well-defined Taylor series expansions about uq. We must now consider two distinct cases: 
the first where \\m.t—,oo{F''^ — G'^) 7^ 0, and the second where this limit vanishes. In the first 
case we must have UxU^ ~ {Fq — GQ)/t^ and so limj^oo^ = 0. At late times we therefore 
expect to recover, at lowest order in 1/t, a comoving spacetime. Physically, this asymptotic time 
evolution just rcficcts simple momentum and angular momentum conservation. An expanding 
area of radius R, velocity U, and mass M oc pR^ will evolve so that MVR is constant; hence 
U (X VUdJ^ oc R-^ oc t-^. 
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If limt^oo(-^'^ — G'"^) = then we have that F' sinhx + G' cosh a:; — *■ 0+ imphes F' cosh a:; + 
G' sinh a; ^ also; hence we must have hmj^oo F' = hm^^oo G" = and so hnif^oo Y = 
tanh(a; + 0{y)) where tanh(^) = limt-KX) G' / F' = G" {uo,y)/F" (uo,y); this wiU be solvable for 
9 because of the requirement that F'^ > G'^ for all finite t. In this second case we should 
therefore expect to recover an asymptotically comoving spacetime only approximately in some 
region about x + 0{y) = 0. 

Consider spacetimes where \imt-^oo{F''^ — G'"^) ^ 0; we find: 



Vo{x,y)Xo{x,y) V tXl Xot 

where the subscript means that a quantity is evaluated at u = uo{x, y), and where: 

r]o{x,y) = Fo{x,y) sinh X + Goix,y) cosh X, 

Xo{x,y) = Fq {x,y) sinh X + GQ{x,y) cosh X, 

vq {x,y) = Fq {x, y) sinh x + Gq (y) cosh x, 

/ ^1 ^ FQ{x,y)F^{x,y) - Go{x,y)G'o{x,y) 
"'^y' F^Hx,y)-G'i{x,y) 

We can see that for such a spacetime to tend to homogeneity at late time, we need 

^ F(,\x,y)-G'i{x,y) 



r]o{x,y)Xoix,y) 



0, 



where a stands for x or y. We shall give an example of a class of spacetimes where this condition 
holds below. 

In the second of the two cases Fq^ — Gq = 0, and asymptotically we then find that either 



cosh 9 (y) 
t'^sinh-'ix + 9{y))F^r,o{x,y) 



if rio{x,y) ^ 0, where rjo is as defined above, or otherwise 

cosh 0(7/) A. 

Kp . ^-^ + 0{t~^), 

t^{-uo{x, y)) sinh-'(x + 9{y))F(lr]o{x, y) 

if r/o = and uq 7^ 0. If 770 = ^^0 = then to leading order the energy density is not positive at 
late times. Wc sec that this class of solutions does not have an isotropic and homogeneous FRW 
limit. We can also see that to leading order it seems that we cannot avoid having a timelike 
singularity at x = —9{y) at late times in this class of solutions. The early-time behaviour of 
both classes of solutions depends strongly on the choice of the free functions F{u, y) and G{u, y). 
Wc shall consider two examples which clearly illustrate the different extremes of behaviour that 
are possible. 



24 



7.3.1 Example 1 

Let us choose F{u, y) = H{y) sinh(u — s{y)) and G{u, y) = H{y) cosh('u — s{y)). We have that: 

J = H{y)smh{x + s{y) -u), 
P = cos\i{x + u — s{y)) — usvahix + u — s{y)) 
= J (^Vl + HH^ - sinh-^(l/iJi) +X + , 



cosh(a; + s{y) — u)P{t, x, y) 

H\y) 
t^P{t,x,y) 

By the analysis of the previous section, this spacctimc has Krolak strong singularities at a; = 
Xo{y, t), where xo(y, t) = sinh^"*^ — \/l + H'^f? — s{y). We note that if i7 > then Kp > for 
all X > xo(y,t)] if H < then Kp > in the region x < xo{y,t); in either case we must restrict 
ourselves only to the region where np is positive. Since P vanishes at x = xq, the circumference 
of the line defined hy x = xq, t = vanishes and this should be properly considered to be a 
single point in the x — y plane. This is a shell-focusing singularity. We now consider the late 
and early-time behaviour of this spacetime. 

Early-time behaviour 

'Early' time now means t « 1/H, and we see that 

P{x, y, t)^t-^{l + In Ht/2 + x + s{y) + 0{{Htf)) 

and 

We can see from this expression that we can see that we will not be able to reach the point i = 
since this leading order term will be become singular at t = to(a^) y) > where: 

to = {2/H) exp (-X - s{y) - 1) > 0. 

We should therefore interpret t = tQ{x,y) as being the true initial singularity - interestingly, 
we note that depending on our choice of s{y), the slice t = to{x,y) can be either spacelike or 
timelike. Ast ^ to, p diverges as {t — to)~^, which is weaker than the we would otherwise expect 
had the fluid not been rotating (i.e. as — ^o)^)- This singularity is of P = type, and so is 
Krolak strong and Tipler weak. 

Late-time behaviour 

'Late' time is H'^t^ S> 1. We can see that 

P{x, y, t) ~ 1 + (x + s{y))/Ht - \{Ht)-^ + 0{{Ht)-^) 
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and: 

At late times (for fixed x and y), this subclass of spacetimes tends to a FRW limit: 
7.3.2 Example 2 

A second class of illustrative spacetimes can we found by taking F{u,y) = C{u) cosh.{9{y)) 
and G{u,y) = C{u) smh.{9{y)). At late times, C'{u) and so these solutions fall into the 
Fq^ — Gq = class mentioned above. The reciprocal time is given by 

^ = C'(n)sinh(x + %)), 

and 

P = (C(u) - uC'{u)) sinh(x + e{y)) 

It is also straightforward to check that the expansion scalar vanishes for this example, 6 = 0, 
and so we must have the shear equal to the vorticity: cr^ = uP" . Let us take, as an example, 

C{u) = ^0 + (C'o — ^o) cosh(M — Mo)> so sinh(u — no) = cosec(x + 6{y))/{Co — AQ)t. We define 
the quantity V := (Co — ^o)^sinh(x + 9) and determine the early and late-time asymptotic 
behaviours. 

Early-time behaviour 

At early times, that is |F| <C 1 for fixed x and y, we have 

u~ - ln|y/2| +uo + 0(F2) 

and 

P^{Ao + i) sinh(x + 9{y)) + ^ In \V/2\ + (1 - uo)/t + 0{V). 
At early times, we there find that the energy density behaves as: 

= ,2 ■ U2r\nf ^^ ii^^ ' ' ^o) + ^(2^0 + 1) sinh(x + 0(y)) + 0{V'))]-' . 

t"^ smh [X + 9{y)) 

Similarly to the previous example is seems as it we will not reach an the point t = 0, as we will 
encountered a singularity at t = to > where ^o > is the value of t for which the quantity 
inside [..] in the above equation vanishes. This singularity is of P = type and so is Krolak 
stronger and Tipler weak. As t ^ to we will again have p oc (t — to)~^. As in example one the 
effect of rotation has been to weaken the strength of the initial singularity. It is important to 
note that t = to(x,y) is not necessarily spacelike. In section 7.5.2 we shall describe an example 
where there is no initial singularity but only a timelike P = type 'central' singularity. 
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Late-time behaviour 

At late times, |y| > 1 we have: n ~ + V'^ + 0{l/V'^) and C(u) ^Cq + F"^(Co - Aq)/2 + 
O^i/V^). Thus, we have 

P - Co sinh(a; + e{y)) - uo/t - {l/2Vt) + 0{1/V^). 

We find for the energy density, assuming Co 7^ that 

~ 7177^ A J ■ >.3, ^., (1 + ^o/(Co sinh(x + e)t)) + 0{t-'). 

t^[Co — ^ojCo smh (a; + U) 

In the case Co = 0, we have 

7.4 A non-vanishing cosmological constant 

Upon the inclusion of a cosmological constant term into the case with one non-zero, non- 
comoving velocity in the x-direction, we find = C{x)sh.p^(t) + D{x)di\{t) where shA(t) = 
A"-^/^ sinh(A^/^t) and chA(t) = cosh(A^/^t). As in the A = case, we can set both Ex and Ey 
to zero. Wc shall define our variables as before: X = K/Jcosh^ 9,Y = tanh^, By = L{x, y, ) and 
E^ = M(x,y,t), Z = XM. With these definitions the equations become: 

M-AM = -ZY'^ 

L = ZY, 

dxM = -L{C{x)sht^{t) + D{x)c\if,{t)), 

Z(C(x)shA(i) + L>(x)chA(i)) = dxL + (C(x)chA(t) + AL>(x)shA(t)) M 

- (C(x)shA(t) + D(x)chA(i))AM. 

When A 7^ 0, we cannot always set = whenever C 7^ without loss of generality; indeed it 
should be noted that those cases where \D\ > |C|/\/A are qualitatively different where the oppo- 
site is true. The former case will generically exhibit a bounce rather than an initial singularity. 
If D^x = 0,, however, then we can, without loss of generality, redefine our x and t coordinates so 
that C = 1 and = 0. We now make the redefinitions T = shA(t)/chA(t) = A"^/^ tanh(A^/^t), 
M = M/chA(t), Z = ZchA(i), and Y = ychA(t). The system {M, L, Z, Y; T} then satisfies the 
A = equations for {M, L, Z, Y; T}. 

7.4.1 Generalisation of the D^x = solution to A / 

We found the general solution for the D = and A = case above. We have just seen that such 
solutions can be easily transformed into A 7^ solutions. We have 

M = shA{t)P{u,x,y), 

where, as before, 

P{u,x,y) = —V? I I — I sinhx-l- I — ) coshx 



u \ u 
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The definition of u in terms of {t, x, y} changes to: 

1 chA(t) 



T shA(t) 



F'(u,y) sinh x + G'{u,y) cosh x . 



The energy density is given hy Kp = Z(l — Y'^)/M. Using our transformation, we have up = 
Z{1 - chX^(t)y)/chi(t)M, and so 

F'-" - + A 

'^^ ch^(t)r2 (F2(t//F)' sinh a; + G2(«/G)' cosh x) (F" sinh x + G" cosh x)' ^ ' 

IF' sinh x + G' cosh x , 

Y = . (114) 

chA {t) F' cosh x + G' sinh x ^ ' 

(115) 

At late times T —>■ and by applying the A = asymptotic analysis we can see that all 

solutions at late times have energy densities that die off in time as e"^^*, and in all cases 
Y ^ as t ^ oo. Thus, at late times all such solutions become comoving iff A > 0. 

7.5 New coordinates for A = case and nature of t = 

When A = the metric of spacetime with C = 1, D = is 

ds^ = -dt^ + t^dx^ + t^P'^dy'^. 

The solutions are equivalent to all C = 0, D = 1, solutions under a coordinate transform: 

t T = t cosh X, x —> X = t sinh X. 

With these new coordinates it is easier to analyse what occurs when t = 0. The metric in {T, X} 
coordinates is: 

ds^ = -dT^ + dX^ + M^dy^ 

where 

M = F{u,y)X + G{u,y)T-u{T,X,y), 
1 = F'{u,y)X + G'{u,y)T. 

The energy density is given by 

p/2 _ QI2 



Kp 



[FX + GT- u){F"X + G"T) ' 



The line i = corresponds to T = ±X, where i — > 0, x ^ xq, xq finite, is T = X = 0, 
X/Y = tanhxQ. This is not necessarily a singularity. 

Consider the behaviour near a point where T = Tq = ±X and Tq finite; wc have I/Tq = 
±F'(uo5 2/) + G'(no, y) and so up < co provided ±F(uo, y) + G(ti+0, y) / uq/Tq and ±F" {uo,y) + 
G"{uo,y) 7^ 0. These conditions together simply require ±F{u,y)+G{u,y) 7^ u/To+o{{u — uo)'^) 
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as u ^ uq. We also require — oo < =fF'(uo, y) + G'{uq, y) < oo. For most choices of F{u, y) and 
G{u, y) the hne T = zizX ^ will, therefore, be non-singular. As we approaches T = X = 0, we 
must have that G'^{u,y) — > oo, and since F'^ > G'^ we conclude that F'"^ must also blow up at 
least as quickly. For many choices of F{u,y) and G{u,y), we fill find that this corresponds to 
u ^ ±oo, and so before we reach the point T = X = OaP = type singularity is encountered, 
where u = FX + GT. 

For now we assume that we can reach the point X = T = by moving along some non- 
spacelike geodesic. We assume that &t X = T = 0, u = uq and that uq is finite. Near u = uowe 
assume that the blow up in F' and G' is due to a pole where F{u, y) A + G{u — uo)~"* and 
G{u, y) ^ B + D{u — uqY, with m > n > — 1, m, n 7^ 0, and m^C^ > v?D'^ ^ 0. So, we have 

nxTiW) • 

where: 

T(X,T,u) = {CX + D{u-uqT-'^T + {AX + BT-uq){u-uq)'^) 
•(m(m + l)CX + n(n + l)D{u - uo)""""?^) 

If m > then the energy density behaves as: 



Kp 



{GX + D{u - uo)"'-"r)(m(m + 1)GX + n{n + l)D{u - uo)"*-"r) 



which clearly blows up as T, X — 0. If — 1 < m < then Kp (x {u — uq)'^X~^, which again 
is manifestly singular. Thus T = X = is, in general, a curvature singularity. Similarly, if 

u ±00 as one approaches T = X = Q and F{u, y) ^ A + Du^, G ^ A + Du^, m > n > 1, 
m?G'^ > ■n?D'^ / 0, then 

m^C^ - n2D2u2(n-m) 

Kp ' 



{GX + Du^-^T){m{m - l)GX + n{n - l)Du'^-^T) ' 

which is manifestly singular at X = T = 0. It is evident from these asymptotics that the LFC 
applies along RNGs terminating at the singularity, but the SFC does not apply. 

As previously stated, it is often the case that the point T = X = is not reachable since it 
lies behind a P = singularity. To illustrate this, and to understand better the nature of these 
spacetimes we construct the Penrose diagrams for two specific choices of F(u,y) and G{u,y). 

7.5.1 Example 1 

We choose F{u, y) = H{y) sinh(u — s(y)), G{u, y) = —H{y) cosh(u — s{y)). This gives: 
u = s{y) + HH{y){X + T)) - ln(l + ^1 + H^y)iT^ - X^)), 

and 

Kp = H\y) [( Vl + i/2(y)(r2 - X2) + u{T, X, y)) ^1 + H^{y){T^ - X^) 
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r=+cio 




Figure 1: Penrose Diagram for the spacetime in Example 1. We have take y fixed and suppressed 
the y-axis. The singuarities are at P = and P' = 0. The fine t = is also drawn. 

There is a spacelike P = 0, Krolak-strong and Tipler-weak singularity when 

siy) + HH{y){X + T)) - ln(l + ^1 + H^iy){T^ - X^)) + ^1 + H^{y){T^ - X^) = 0, 

and a timelike P' = 0, Krolak weak sing ularity at X^ = 1/H^{y) + T^. The point X = T = Q lies 
beyond the boundary of this spacetime. This space is homogeneous at late times: Kp ~ l/T^, 
for fixed y and X. We construct the Penrose diagram of this spacetime for fixed y in figure 
1. We can see from fig. 1 that all past-directed timelike geodesies terminate on a singularity, 
either on the 'big-bang', singularity at P = 0, or on the timelike P' = line. Although we 
have labelled the P = singularity as a 'big-bang' it is important to note that it need not be 
everywhere spacclikc, and different choices of H{y) and s(y) can easily result in it being timelike 
in some locality - that is appears spacclikc in fig. 1 is due to our suppression of the y-axis. Since 
the singularity at P' = is weak in the sense of Krolak it is, in principle, possible to continue 
through it, however we shall not consider here what may lie beyond it. 

7.5.2 Example 2 

We take F = C(y)e"cosh(0(y)), G = C(y)e" sinh(0(y)), so 

u = -\n{C{y)Xcosh{e{y)) + C(y)rsinh(^(y))). 

In this case 

(l-n) 
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r=+oo 



r 




r 

r=-oo 



Figure 2: Penrose Diagram for the spacetime in Example 2. We have take y fixed and suppressed 
the y-axis. The only singularity is of P = type. The lines t = and T = are also drawn. 

and the only singularities are of P = type and occur when u = 1. Unlike in the previous 
example, the P = singularity is timelike in this case, and can be thought of as a centre. At 
late times the space is not homogeneous and Kp on 1/ (T^ In T) . As before, the point T = X = 
lies beyond the boundary of the spacetime. We construct the Penrose diagram for this spacetime 
(at fixed y) in figure 2, from this it is clear that there is no 'big-bang' initial singularity in this 
model. Indeed it can be easily checked that the expansion scalar, 0, vanishes in for this solution 
and so this is actually an example of an inhomogeneous static spacetime with non-vanishing 
rotation and shear. Using the results of sections 3.2 3.3 we see that we must have = uP' and 
71 = 2Kp where: 



8 Scalar-Field Spacetimes 

8.1 Solutions with one spacelike Killing vector 

We now find the general solution to 2+1 gravity with a massless scalar- field source under the 
assumption that there is one spacelike Killing vector. The metric takes the following form: 

ds^ = -2A{u, v)dudv + C^(tt, v)dy'^. 
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We can transform this to t, x coordinates by defining u = {t+x)/ v^, v = (t — x)/ \pl. We assume 
a scalar field source, = c(>{u,v), witfi energy momentum tensor: Tab = da4>df,(^ — ^gab{d4')^ 
where (p satisfies the conservation equation 

□0 = 0. 

With these prescriptions the (yu) and {yv) components of Einstein equations are satisfied 
trivially, and the {yy) component requires C^uv = Oj the general solution of which is: 

C{u,v) = f{v)+g{u) 

We also define D{u,v) = f{v) — g{u), and move from {u,v) to (C, -D) coordinates. By making 
different choices of f{v) and g{u) we can arrange that either: daC is spacelike and daD is 
timclikc, or daD is spacelike and daC is time-like, or both daC and daD are null. In terms of C 
and D the cp wave equation reads: 

(l>,DD = ^,CC + 4>,c/C. 

This is just the wave-equation in cylindrical polar coordinates with axial and azimuthal symme- 
try (with D playing the role of the usual time coordinate and C of the radial coordinate). We 
can solve this in terms of Bessel functions. 




where A{k), B{k), E{k) G C are arbitrary and Jq and Yq are zero-order Bessel functions of the 
first and second kind respectively. Finally, we solve the equations for A to give: 

1^ ( fn = 2^ r ^D'M,D + FiC) 
\f'(v)g'{u)J Jd^^c) 

where ^ 

F{C)= [ dC'C'{^% + ^%)\n=Do(C'y 

J Co 

Boundary conditions for (f) are need to specify the solution further. 
8.2 PP-wave spacetimes 

In 2-1-1 spacetimes the metric for a scalar-field PP-wave spacetime can be written in the form: 

ds^ = H{u, x)du'^ + 2dudv + da;^. 

The Einstein equations read Ruu = — l/2-ff,a;a: = n{4',u)'^ and Rab = otherwise; this implies 
(f) = (f>{u), where (l){u) is arbitrary. Solving for H, we find: 

H = A{u) + B{u)x - K{^'iu)f x^, 

where A{u) and B{u) can both be freely specified. It can be checked that in 2+1 dimensions 
these are the only perfect-fiuid solutions (they are equivalent to an irrotational p = p fiuid) 
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that are compatible with the PP-wave metric ansatz given above. Indeed, we find that the only 
permitted choices of energy-momentum tensor must satisfy T^u 7^ HT and Tab = Tgab otherwise. 
We find similar PP-wave solutions by considering the 2+1 Einstein-Maxwell equations. Up to 
gauge transformations, all the solutions are of the form Au = Ay = 0, and A^ = 4>{u) for 
the electromagnetic potential, with (l){u) arbitrary and H{u, x) as given above. The only non- 
vanishing components of F^b are F^^ = —Fxu = 4',u{u). 

9 Discussion 

Whereas the general cosmological solutions of the 3-1-1 dimensional Einstein equations are in- 
tractably complicated and likely dominated by non-integrability, the structure of the theory in 
2+1 offers the possibility of making considerable progress towards finding the general solution 
in several interesting situations. This fact, together with our current perception that quantum 
field theory fits more naturally in three rather than four dimensions, has motivated the study 
of Einstein's theory in 3-dimensional spacetimes. 

In this article we employed covariant and first-order formalism techniques to study the 
properties of general relativity in three dimensions. The covariant approach provided an ir- 
reducible decomposition of the relativistic equations and allowed for a mathematically compact 
and physically transparent description of their properties. Using this information we reviewed 
the kinematical, dynamical and geometrical features of 3-dimensional spacetimes and identified 
the special features that distinguish them from the standard 3+1 models. These include the 
key role of the isotropic pressure as the sole contributor to the gravitational mass of the system 
and the fact that vorticity never increases with time. We also reviewed the 3-d analogues of 
the spatially homogeneous and isotropic FRW models and investigated their stability against 
linear perturbations. We found that, unlike their conventional counterparts, dust-dominated 3-d 
homogeneous and isotropic spacetimes are stable under shear and vorticity distortions and (neu- 
trally) stable against disturbances in their density distribution. The latter reflects the vanishing 
of the total gravitational mass in 3-dimensional dust models, which ensures the absence of linear 
Jeans-type instabilities. In addition to isotropic spacetimes, we also looked at 3-dimensional 
anisotropic models providing Kasner-like solutions for the case of pressure-free matter and gen- 
eralising Godcl's universe to three dimensions. The covariant formalism allowed us to carry out 
these analyses by a study of the kinematic variables characterising the expansion of the universe. 
The absence of both electric and magnetic Weyl curvature components in three dimensions con- 
siderably simplifies the analysis. We then specialised further to the case of a pressureless matter 
source. In addition to being physically realistic, this assumption produces a significant further 
simplification of the cosmological field equations in three-dimensional spacetimes. We were able 
to find the general cosmological solutions of the theory in the case where the matter was comov- 
ing. No symmetry assumptions were made. We then considered the fully general pressureless 
fluid system with non-comoving velocities. We were able to solve the system in the case where 
one spatial velocity component was zero whilst the other was non-zero. This allowed us to carry 
out an asymptotic study, close to and far from singularities, of an inhomogeneous cosmology 
with rotation, expansion and shear. All the singularities arising in these solutions were classified 
using the different criteria of strength introduced by Krolak and Tipler. We were able to provide 
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a simple transformation which generahsed all the solutions we found with vanishing cosmological 
constant into new solutions with non-zero cosmological constant. Finally, we considered scalar- 
field metric with one Killing vector and found all the PP-wave solutions in 2+1 dimensional 
universes. 

These investigations suggestion a number of problems for further study. Exact solutions in 
the cases with non-zero isotropic and anisotropic pressure remain to be investigated. In the case 
of zero pressure, we have analysed the problem of the general solution of the three-dimensional 
Einstein equations into a well-defined system of partial differential equations. We have solved for 
the case with comoving velocities and a single non-comoving velocity but the problem remains 
to find the general solution of the equations when both non-comoving fluid velocities are present. 
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